Методы решения дифференциальных уравнений
1. Уравнения с разделяющимися переменными
\[
\frac{dy}{dx} = f(x)g(y)
\]
Алгоритм решения:
1. Записать уравнение в виде: \( \frac{dy}{dx} = f(x)g(y) \)
2. Разделить переменные: \( \frac{dy}{g(y)} = f(x)dx \)
3. Проинтегрировать обе части: \( \int \frac{dy}{g(y)} = \int f(x)dx + C \)
4. Выразить y через x (если возможно)
Пример:
Решить уравнение: \( y' = x y^2 \)
1. Записываем в виде: \( \frac{dy}{dx} = x y^2 \)
2. Разделяем переменные: \( \frac{dy}{y^2} = x dx \)
3. Интегрируем: \( \int y^{-2} dy = \int x dx \)
4. Получаем: \( -\frac{1}{y} = \frac{x^2}{2} + C \)
5. Выражаем y: \( y = -\frac{1}{\frac{x^2}{2} + C} = -\frac{2}{x^2 + 2C} \)
Особые случаи:
- Если \( g(y) = 0 \) при некотором \( y = y_0 \), то \( y = y_0 \) - особое решение
- Если \( f(x) \) или \( g(y) \) не интегрируются в элементарных функциях, решение может быть выражено в неявной форме
2. Однородные уравнения
\[
y' = f\left(\frac{y}{x}\right)
\]
Алгоритм решения:
1. Записать уравнение в виде: \( y' = f\left(\frac{y}{x}\right) \)
2. Сделать замену: \( u = \frac{y}{x} \), тогда \( y = ux \), \( y' = u'x + u \)
3. Подставить в уравнение: \( u'x + u = f(u) \)
4. Разделить переменные: \( \frac{du}{f(u) - u} = \frac{dx}{x} \)
5. Проинтегрировать: \( \int \frac{du}{f(u) - u} = \int \frac{dx}{x} + C \)
6. Вернуться к исходным переменным, подставив \( u = \frac{y}{x} \)
Пример:
Решить уравнение: \( y' = \frac{y}{x} + \tan\left(\frac{y}{x}\right) \)
1. Замена: \( u = \frac{y}{x} \), \( y = ux \), \( y' = u'x + u \)
2. Подставляем: \( u'x + u = u + \tan(u) \)
3. Упрощаем: \( u'x = \tan(u) \)
4. Разделяем переменные: \( \frac{du}{\tan(u)} = \frac{dx}{x} \)
5. Интегрируем: \( \int \cot(u) du = \int \frac{dx}{x} \)
6. Получаем: \( \ln|\sin(u)| = \ln|x| + C \)
7. Возвращаемся к y: \( \sin\left(\frac{y}{x}\right) = Cx \)
Как определить однородное уравнение:
Уравнение \( P(x,y)dx + Q(x,y)dy = 0 \) является однородным, если функции P и Q однородны одной степени:
\( P(tx,ty) = t^n P(x,y) \), \( Q(tx,ty) = t^n Q(x,y) \)
3. Линейные уравнения
\[
y' + P(x)y = Q(x)
\]
Алгоритм решения:
1. Записать уравнение в стандартной форме: \( y' + P(x)y = Q(x) \)
2. Найти интегрирующий множитель: \( \mu(x) = e^{\int P(x)dx} \)
3. Умножить уравнение на интегрирующий множитель:
\[
\mu(x)y' + \mu(x)P(x)y = \mu(x)Q(x)
\]
4. Левая часть является производной: \( \frac{d}{dx}[\mu(x)y] = \mu(x)Q(x) \)
5. Проинтегрировать: \( \mu(x)y = \int \mu(x)Q(x)dx + C \)
6. Выразить y: \( y = \frac{1}{\mu(x)} \left( \int \mu(x)Q(x)dx + C \right) \)
Пример:
Решить уравнение: \( y' + 2xy = x \)
1. Стандартная форма: \( y' + 2xy = x \) (уже в нужной форме)
2. Интегрирующий множитель: \( \mu(x) = e^{\int 2x dx} = e^{x^2} \)
3. Умножаем: \( e^{x^2}y' + 2xe^{x^2}y = xe^{x^2} \)
4. Левая часть: \( \frac{d}{dx}[e^{x^2}y] = xe^{x^2} \)
5. Интегрируем: \( e^{x^2}y = \int xe^{x^2} dx = \frac{1}{2}e^{x^2} + C \)
6. Решение: \( y = \frac{1}{2} + Ce^{-x^2} \)
Особые случаи:
- Если \( Q(x) = 0 \), уравнение называется однородным линейным уравнением
- Решение однородного уравнения: \( y = Ce^{-\int P(x)dx} \)
4. Уравнения Бернулли
\[
y' + P(x)y = Q(x)y^n
\]
Алгоритм решения:
1. Записать уравнение в стандартной форме: \( y' + P(x)y = Q(x)y^n \)
2. Разделить на \( y^n \) (при \( y \neq 0 \)): \( y^{-n}y' + P(x)y^{1-n} = Q(x) \)
3. Сделать замену: \( z = y^{1-n} \), тогда \( z' = (1-n)y^{-n}y' \)
4. Подставить в уравнение: \( \frac{z'}{1-n} + P(x)z = Q(x) \)
5. Умножить на \( (1-n) \): \( z' + (1-n)P(x)z = (1-n)Q(x) \)
6. Решить как линейное уравнение относительно z
7. Вернуться к переменной y: \( y = z^{\frac{1}{1-n}} \)
Пример:
Решить уравнение: \( y' + \frac{y}{x} = x^2 y^3 \)
1. Стандартная форма: \( y' + \frac{1}{x}y = x^2 y^3 \)
2. Делим на \( y^3 \): \( y^{-3}y' + \frac{1}{x}y^{-2} = x^2 \)
3. Замена: \( z = y^{-2} \), \( z' = -2y^{-3}y' \)
4. Подставляем: \( -\frac{1}{2}z' + \frac{1}{x}z = x^2 \)
5. Умножаем на -2: \( z' - \frac{2}{x}z = -2x^2 \)
6. Решаем линейное уравнение:
Интегрирующий множитель: \( \mu(x) = e^{\int -\frac{2}{x} dx} = e^{-2\ln|x|} = \frac{1}{x^2} \)
\( \frac{1}{x^2}z' - \frac{2}{x^3}z = -2 \)
\( \frac{d}{dx}\left[\frac{z}{x^2}\right] = -2 \)
\( \frac{z}{x^2} = -2x + C \)
\( z = -2x^3 + Cx^2 \)
7. Возвращаемся к y: \( y^{-2} = -2x^3 + Cx^2 \)
8. Решение: \( y = \frac{1}{\sqrt{-2x^3 + Cx^2}} \)
Особые случаи:
- При n = 0 уравнение становится линейным
- При n = 1 уравнение становится линейным (делим на y)
- Решение y = 0 всегда является особым решением при n > 0
5. Уравнения в полных дифференциалах
\[
M(x,y)dx + N(x,y)dy = 0
\]
где \( \frac{\partial M}{\partial y} = \frac{\partial N}{\partial x} \)
Алгоритм решения:
1. Проверить условие: \( \frac{\partial M}{\partial y} = \frac{\partial N}{\partial x} \)
2. Найти функцию U(x,y) такую, что:
\[
\frac{\partial U}{\partial x} = M(x,y), \quad \frac{\partial U}{\partial y} = N(x,y)
\]
3. Проинтегрировать M по x: \( U(x,y) = \int M(x,y) dx + \varphi(y) \)
4. Найти \( \varphi(y) \) из условия: \( \frac{\partial U}{\partial y} = N(x,y) \)
5. Общее решение: \( U(x,y) = C \)
Пример:
Решить уравнение: \( (2xy + y^2)dx + (x^2 + 2xy)dy = 0 \)
1. Проверяем условие:
\[
M = 2xy + y^2, \quad N = x^2 + 2xy
\]
\[
\frac{\partial M}{\partial y} = 2x + 2y, \quad \frac{\partial N}{\partial x} = 2x + 2y
\]
Условие выполняется.
2. Находим U(x,y):
\[
U(x,y) = \int (2xy + y^2) dx = x^2y + xy^2 + \varphi(y)
\]
3. Дифференцируем по y:
\[
\frac{\partial U}{\partial y} = x^2 + 2xy + \varphi'(y)
\]
4. Приравниваем к N:
\[
x^2 + 2xy + \varphi'(y) = x^2 + 2xy \Rightarrow \varphi'(y) = 0 \Rightarrow \varphi(y) = C_1
\]
5. Общее решение: \( x^2y + xy^2 = C \)
Если уравнение не является уравнением в полных дифференциалах:
Можно попытаться найти интегрирующий множитель μ(x,y) такой, что:
\( \frac{\partial (\mu M)}{\partial y} = \frac{\partial (\mu N)}{\partial x} \)
Шпаргалка по методам решения
Определение типа уравнения:
- Попытаться разделить переменные
- Проверить на однородность
- Проверить, является ли уравнение линейным
- Проверить, является ли уравнение уравнением Бернулли
- Проверить, является ли уравнение уравнением в полных дифференциалах
Полезные замены:
- \( u = \frac{y}{x} \) для однородных уравнений
- \( z = y^{1-n} \) для уравнений Бернулли
- \( y = ux \) для приведения к однородному виду
Часто встречающиеся интегралы:
\[
\int \frac{dx}{x} = \ln|x| + C
\]
\[
\int e^{ax} dx = \frac{1}{a} e^{ax} + C
\]
\[
\int \sin(ax) dx = -\frac{1}{a} \cos(ax) + C
\]
\[
\int \cos(ax) dx = \frac{1}{a} \sin(ax) + C
\]
\[
\int \frac{dx}{x^2 + a^2} = \frac{1}{a} \arctan\left(\frac{x}{a}\right) + C
\]
Методы решения уравнений высших порядков:
- Понижение порядка для уравнений вида \( y^{(n)} = f(x) \)
- Замена \( y' = p \) для уравнений, не содержащих y явно
- Замена \( y' = p \), \( y'' = p \frac{dp}{dy} \) для уравнений, не содержащих x явно
- Линейные однородные уравнения с постоянными коэффициентами
- Метод вариации постоянных для неоднородных уравнений